Solving three-body scattering problem in the momentum lattice representation 
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A brief description of the novel approach towards solving few-body scattering problems in a finite- 
dimensional functional space of the Z/2-type is presented. The method is based on the complete few- 
body continuum discretization in the basis of stationary wave packets. This basis, being transformed 
to the momentum representation, leads to the cell-lattice- like discretization of the momentum space. 
So the initial scattering problem can be formulated on the multi-dimensional momentum lattice 
which makes it possible to reduce the solution of any scattering problem above the breakup threshold 
(where the integral kernels include, in general, some complicated moving singularities) to convenient 
simple matrix equations which can be solved on the real energy axis. The phase shifts and inelasticity 
parameters for the three-body nd elastic scattering with MT I-III NN potential both below and 
above the three-body breakup threshold calculated with the proposed wave-packet technique are in 
a very good agreement with the previous accurate benchmark calculation results. 
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1. Motivation of the work . The strictly proved inte- 
gral equations for the solution of few-body scattering 
problems were developed many years ago by Faddeev and 
Yakubosky After these pioneer works a lot of in- 

vestigations in the few-body quantum physics have been 
done along these lines for a few last decades. In spite of a 
great progress in this field [1, 0, IE IE 0, II, Q , the practical 
solution of the few-nucleon scattering problems with real- 
istic 2N- and 3N- interactions, especially above the three 
body breakup threshold, is still remained rather cum- 
bersome computational problem which needs an appeal 
to very powerful computer resources. Moreover, even 
till now AN systems above the three-body threshold can 
be practically treated within the Faddeev- Yakubovsky 
framework only with simple pairwise local interactions 
0, Q . The reason is in very laborious numerical routines 
in the coordinate space or complicated moving singulari- 
ties in kernels of the integral momentum-space equations. 

At the same time, several efficient methods for the 
approximation of few-body continuum wavefunctions in 
various Lo bases have been developed [l(| EH 
H EH ElES El- These are the "moment T-matrix 
method" [lfj|. J-matrix approach [ill . H3. "the har- 
monic oscillator representation" [14j |. the Lorentz inte- 
gral transform method [l5j . the continuum-discretized 
coupled-channel method (CDCC) [M El El etc. How- 
ever most of them can be used for special cases of the 
few-body scattering only, e.g. for the so called truly few- 
body scattering when there are no bound states in any 
two-body subsystems [lj], or for the composite particle 
scattering off heavy target when stripping channels can 
be neglected 0, ES EH- In other cases one describes 
the processes when the few-body wavefunctions in the 
initial channel are of bound-state type and L2 basis is 



used to approximate the final-state few-body continuum 
only [l5j], or one treats a three-body scattering at small 
energies below the three-body threshold only [l3j]. So 
that, with the above L2-type methods no precise calcula- 
tions for the basic three-body n — d scattering case above 
the breakup threshold have been carried out up to date. 1 

Thus, it would be very convenient to have in our dis- 
posal a sufficiently universal method for general contin- 
uum discretization in different two- and few-body scatter- 
ing problems (in nuclear, atomic, hadronic etc. physics), 
which operating with L2 functions only and non-singular 
matrix equations both below and above the breakup 
thresholds. 

Few years ago the present authors have developed a 
new approach to solving few-body scattering problems 
based on discretization of continuous spectrum of total 
Hamiltonian [H EH, [H [H, [I| . The method uses the 
stationary wave packets, i.e. are L2-type functions, in- 
stead of the exact scattering wave functions. In these 
works an original wave-packet formalism have been de- 
veloped which allowed to construct finite-dimensional (f.- 
d.) approximations for basic scattering-theory operators 
and find the scattering observables using such approx- 
imations. The approach have recently been tested for 
the elastic scattering and breakup of a composite projec- 
tiles scattered off heavy targets (with neglecting the strip- 
ping processes), and a perfect agreement with the CDCC 
method results have been found [23, [24( . In the present 
paper we extend our wave-packet approach much further, 
towards solving a general three-body scattering problem 
on the base of the projected Faddeev equations and illus- 
trate it on example of quartet and doublet n— d scattering 
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The realistic three-nucleon calculations for e.g. nd (or pd) scat- 
tering below and above three-body breakup threshold have been 
done either with variational method [OJ using the Schrodinger 
equation approach or with the Faddeev equations in the momen- 
tum [§[ or in the configuration space 0] . 
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below and above three- nucleoli breakup threshold 2 . 

2. Formulation of the approach. Here we describe the 
three-body wave-packet discretization procedure for the 
elastic n — d scattering. The elastic amplitude X for the 
quartet case can be found from a single integral Faddeev 
equation [25j 



X = -Pvt - PviGxX, 



(1) 



where V\ is the triplet NN interaction potential, G\ — 
(E — Hi) -1 is the three-body resolvent of the channel 
Hamiltonian Hi = Ho + v\ and P is the permutation 
operator. 

Let us introduce some finite basis {\Si)} J ^ 1 such that 
the projector T^f onto the f.-d. basis subspace can 
(in some not rigorous sense) approximate the unit op- 
erator rV = YliLi \Si)(Si\ — > 1. Then one can de- 
fine the A/'-dimensional approximation for any operator 
A as its projection onto the respective basis subspace: 
T^fATj^f = Yljj \Si)Aij(Sj\, with corresponding matrix 
elements Ay = (Si\A\Sj). Using such matrix approxi- 
mations for the scattering operators, the initial integral 
equation can be reduced to the respective matrix equa- 
tion. We will denote the matrices of projected operators 
with corresponding bold letters. So, one finds the matrix 
equation instead of the integral equation ([1]): 



X 



Pvi - PviGiX. 



(2) 



Thus, it looks like it would be possible to find some ap- 
proximate solutions of the initial integral equation us- 
ing some appropriate Li bases through the simple ma- 
trix algebra. However, not everyone Li basis suits for 
this purpose. The integral kernel of the Faddeev equa- 
tion includes the fixed-pole singularities and also the 
complicated moving singularities above the three-body 
breakup threshold. Just these singularities correspond 
to the proper boundary conditions in coordinate space 
and provide the correct physical solution of the Faddeev 
equations, but a construction of the appropriate basis for 
the projection of such kernels is highly non-trivial prob- 
lem. Another key problem here is a calculation of matrix 
elements for v\ , P and especially for G\ operators in the 
chosen basis, that can appear, in general, rather cum- 
bersome task. In this point the appropriate stationary 
wave-packet basis helps to overcome all the above diffi- 
culties and formulate the matrix equations whose solu- 
tions can really approximate the solutions of the initial 
integral equations. Such a basis will be demonstrated to 
provide a convenient analytical diagonal representation 
for the three-body channel resolvent matrix Gi and, on 
the other hand, this basis covers a sufficiently wide func- 
tional Li space to provide well converged results. 



2 As far as the present authors are aware this is the first precise 
fully L2 approximated solution for the Faddeev equation above 
the three-body threshold. 



The channel Hamiltonian Hi is the direct sum of two 
subHamiltonians corresponding to the system motion 
along two independent Jacoby coordinates: 



Hi = hi 8 h , 



(3) 



where subHamiltonian hi defines the NN subsystem mo- 
tion (including interaction vi) and subHamiltonian ho 
corresponds to the free motion of the third nucleon rel- 
atively the NN subsystem center of mass. Now let's 
define two sets of momentum bins [pm-i,Pm}m=i and 
[q n -i, q n ]n=i corresponding to the continuum discretiza- 
tions for subHamiltonians hi and ho- The (two-body) 
stationary wave-packet bases (WPB) are defined as inte- 
grals of exact continuum wave functions \ip p ) and \tpo q ) of 
subHamiltonians hi and ho over the respective momen- 
tum bins: 



f 



dp\ip p ) 



1 



\JJn 



dq\ip, 



'Oqi 



<Zn- 



'dn 

(4) 

where b m = p m ~ Pm-i and d n = qn~1n-i are bin widths. 
Now, the three-body WPB states \Sij) are defined just as 
products of two-body wave-packet states \zi) (including 
the bound state wave function |zq)) and \yj): 



\Sij) = 



(5) 



We omit here partial wave indices for the sake of simplic- 
ity. The properties of the wave-packet sets \zj) and \yj) 
have been investigated in detail in refs. [2(| El, H2, [23 |. 
In particular, the very useful property of such a packet 
basis is that the matrices for the projected resolvents 
of subHamiltonians gi(E) = (E + iO — /ii) 1 and go — 
(E + iO — /»o) 1 are diagonal and defined by a simple an- 
alytical expressions depending only on the spectrum dis- 
cretization parameters. This property of WPB allows us 
to construct the f.-d. analytical diagonal representation 
for the channel resolvent Gi(E), which is a convolution 
of two-body resolvents gi and go- 

Indeed, the exact three-body channel resolvent can be 
written [2lH as a sum of two terms Gi(E) = Gf G (E) + 
Gi°, where the bound-continuum part has the form: 



/>oo 

Gf C (E) = / dq 
Jo 



\Zo,1pOq}(Zo,1pQg\ 



E + iO-e Q 



3q 2 ' 
4m 



(6) 



and eo is the binding energy for the (single) NN bound 
state. The continuum-continuum part takes the form: 



/•oo roo 

G° C (E) = dp dq 
Jo Jo 



_ 3q 2 
4m 



E + iO 



(7) 



Now let's construct a projection of the exact channel re- 
solvent onto the three-body WPB. The following analyt- 
ical formulas for the diagonal f.-d. projection of Gi can 
then be found: 



^ j G^\zo,y j ){zo,yj\, 



G ?°= E GfP\z hVj ){z h yA 



(8) 
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where the matrix elements G^p and Gfj in eq. ([8]) are 
defined as integrals over the respective momentum bins: 



1 



dq 



3q 2 5 
4 m 



G 



cc 



didj 



dpdq 



These matrix elements depend, in general, on the spec- 
trum partition parameters (i.e. pi and qj values). How- 
ever we found when the wave-packet expansions of the 
three-body amplitude is convergent the final result turns 
out to be independent upon the particular spectral par- 
tition parameters. The integrals in eqs.([8K) and (03 ) are 
analytically tractable that gives a simple analytical f.-d. 
representation for the three-body channel resolvent G\. 
Such an analytical representation is main feature of the 
wave-packet approach since it allows to simplify solution 
of the general three-body scattering problem drastically. 

Now the key question arises: how to construct practi- 
cally the above wave-packet basis. The free packets \y n ) 
corresponding to the free motion of third nuclcon rela- 
tively the NN subsystem cm. take in the momentum 
representation the form of simple step-like functions: 



vim) 



%-<?j-x) -6{q-qj) 



(9) 



where 6{x) is the Heaviside function. 

The scattering wave packets \zi) describing the scat- 
tering in the NN two-body subsystem can be very well 
approximated by pseudostates \zi) obtained from the di- 
agonalization of the subHamiltonian hi in some appro- 
priate L2 basis 20]. In the present work we use for this 
diagonalization a free wave-packet basis \xk) correspond- 
ing to the free NN motion. Thus, we solve the two-body 
variational problem directly on the free WPB and as a 
result obtain a set of variational functions 



M 

E 

k=Q 



O lk \x k ) 



0, . . . , M, 



(10) 



the first of which (for the problem in question) is the 
wavefunction of the bound state (deuteron) and other 
ones are very good approximations for the exact scatter- 
ing packets. 

In the Fig. 1 the coordinate-space behavior of some 
from the first 50 variational functions (including the 
deuteron) is shown. It is clear that the free packet ba- 
sis allows to approximate the respective scattering wave 
functions up to a very far asymptotic region (in the Fig- 
ure the functions Zi(r) are given at r < 80 fm, but actu- 
ally they coincide with exact scattering wave packets up 
to r ~ 1000 fm). This long-range behavior of the basis 
functions plays a crucial role in the three-body scatter- 
ing, especially above the breakup threshold, because it 
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FIG. 1: (Color online) Bound state (solid curve) and several 
interaction wave packets (different dashed lines) for the MT 
III potential constructed from the free momentum packets in 
the coordinate space. 



provides a proper overlapping between basis functions in 
different Jacoby-coordinate sets. It should be mentioned 
that in our previous calculations plj we used Gaussian 
basis to approximate the interaction packets in the whole 
space and so the wrong long-range behavior of the basis 
functions did not allow us to obtain well converged re- 
sults above the break-up threshold (while for the smaller 
energies the Gaussian approximation works well) 3 . 

Besides, the momentum lattice basis is very convenient 
to find the matrix elements of the permutation operator 
P. Using approximation (fTOf for the scattering packets 
\zi), these matrix elements can be expressed in the form: 



E 

kk' 



OikO*, k ,P®. 



kj,k'j' 



(11) 



where P k jk'j' = { x kVj\P\xk'yj') is the permutation ma- 
trix elements taken on the two-dimensional free wave 
packets (for a two-dimensional lattice). Using hy- 
perspherical momentum coordinates, the calculation of 
P®j k , y can be reduced to a one-dimensional numerical 
integration over hypermomentum p 2 + q 2 ■ The tech- 
nique of this calculation will be given in detail else- 
where. It should be stressed here that this is one of 
the key points for the whole our approach. In fact, 
to solve two-dimensional Faddeev equations by conven- 
tional methods [H one needs (due to the appearance of 
permutation operator P in the integral kernel) to use a 
very time-consuming two- and three-dimensional interpo- 
lations (many thousands or even millions of such interpo- 



s It should be noted the free packet basis (i.e. the step-like func- 
tions in the momentum space) is, of course, not optimal for cal- 
culation of bound states. E.g. in our case only 20 Gaussian 
functions are necessary to obtain Ei, = —2.225 MeV and ca. 100 
step-like functions are required to reach the same precise bound 
energy value. However very good approximation of scattering 
wave functions in two-body subsystems is the decisive factor here. 
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lations) at each iteration step to find the solution in the 
initial Jacoby set from the "rotated" (by the permutation 
P) Jacoby sets. So, such numerous multi-dimensional 
interpolations at each step take a big portion of compu- 
tational time in practical solutions of three-body inte- 
gral equations. When solving the four-body Yakubovsky 
equations the dimension and number of the above each 
step interpolations gets even higher. Thus, the wave- 
packet approach allows to avoid such multi-dimensional 
interpolations. 

After solving the matrix eq. the on-shell elas- 

tic amplitude A e \(E) in wave-packet approximation can 
be found as a diagonal (on-shell) matrix element of X- 
matrix (which is a solution of the matrix equation |[2}): 



2m X 0no 



On a 



3<?o 



(12) 



where index no being denote the singular q-b'm to which 

the on-shell momentum qo = yj ^m(E — eg) belongs: 

<Zo £ (Sn -i> Qn )- Let's notice that in order to find the 
elastic amplitude according to eq. (fT2Tj , one needs to solve 
one linear equation only just for one column X mn $ nQ1 but 
not to do a full matrix inversion in eq. ([2j). 

3. Numerical results . To illustrate the accuracy and 
effectiveness of the proposed wave-packet technique we 
calculated the real phase shifts and inelasticity parame- 
ters for the three-body elastic n—d scattering in the quar- 
tet and doublet S'-wave channels with the model Malfliet- 
Tjon NN potential MT-III. The results of these calcula- 
tions are shown in Figs. 2 and 3 for the spin-quartet and 
in Figs. 4 and 5 for the spin-doublet channels respectively. 
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FIG. 2: (Color online) The energy dependence of the real 
phase shift for S-wave quartet nd scattering calculated by 
means of the momentum-packet discretized Faddeev equation 
at different dimensions M x N of the lattice basis: 100 x 100 
(dashed curve), 200 x 200 (solid curve). Results of the direct 
Faddeev equation solution from ref.0, @] are marked as 



In the case of doublet scattering one has the system of 
two matrix equations instead of one matrix equation @ 
where two amplitudes correspond to two possible spin 
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FIG. 3: (Color online) The same as in Fig. [2] but for the 
inelasticity parameter n. 



states (triplet and singlet) of the NN subsystem. To 
check the accuracy of the method we have compared our 
results with the previous benchmark calculation results 
from ref. Q (below the deuteron breakup threshold) and 
ref. [|| (above the deuteron breakup threshold) The exact 
parameters of the NN potential are taken from ref. 0]. 
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FIG. 4: (Color online) The energy dependence of the real 
phase shift for the 5-wave doublet nd scattering calculated 
by means of momentum-packet discretized Faddeev equation 
at different dimensions M x N oi the lattice basis: (50 + 
50) x 50 (dashed curve) , (80 + 80) x 80 (dotted curve) and 
(100 + 100) x 100 (solid curve). Results of the direct Faddeev 
equation solution from ref. 0,0] are marked as A. 

As it is seen from the Figs. 2-5 the wave-packet dis- 
cretization technique for the three-body continuum works 
successfully for the general three-body scattering prob- 
lem both below and above breakup threshold. Thus, 
at the first time we have solved the three-body scatter- 
ing problem above the breakup threshold using only f.- 
d. approximation of the L2-type for the Faddeev kernel. 
Just the use of the momentum-lattice basis allowed us to 
achieve a good convergence and accuracy on this way. It 
is interesting to remark that although the quartet case 
seems to be simpler from the first glance (one equation 
and the inelasticity is less than that in the doublet case) 
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it turns out to be more difficult numerical problem (the 
larger basis is needed for convergence) for the wave packet 
discretization approach. 
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FIG. 5: (Color online) The same as in Fig. [4] but for the 
inelasticity parameter r\. 

4. Conclusion. Let's briefly outline here the most im- 
portant points of this study. For the first time the three- 
body scattering problem in the Faddeev framework above 
the breakup threshold has been successfully solved in the 
three-body L2 basis representation using the lattice ap- 
proximation scheme (which is the technique of the three- 
body continuum discretization) . The success and advan- 
tages of the lattice approach are related to the following 
distinctive features. 

(i) The explicit analytical f.-d. approximation for the 
three-body channel resolvent G\ allows to reduce initial 
integral Faddeev equation to the matrix one that can be 



solved directly on the real energy axis. 

(ii) The scattering wave packets (corresponding to the 
NN interaction) can be approximated by pseudostates 
of the two-body NN subHamiltonian matrix in the free 
wave-packet basis, which allows to avoid calculation of 
the two-body i-matrix and obtain explicitly matrix ele- 
ments of the permutation operator P that includes over- 
lapping between wave-packet basis states of the different 
three-body channel Hamiltonians. 

(iii) This convenient closed form for the matrix el- 
ements of the permutation operator P in the WPB 
also makes it possible to avoid completely very time- 
consuming multi-dimensional interpolations of the iter- 
ated kernels which are usually assist in conventional tech- 
niques of the Faddeev equation numerical integration in 
the momentum space. 

(iv) The very long-range type of the wave-packet func- 
tions (non-vanishing at distances ~ 1000 fm) allows 
to approximate properly the overlapping between basis 
states in different Jacoby coordinate sets. This leads to 
the proper asymptotic behavior of the solutions along dif- 
ferent Jacoby coordinates, which couldn't be provided by 
means of conventionally used short-range type L2 bases. 

Besides, this long-range behavior of the wave-packet 
basis functions looks also very promising for the proper 
incorporation of the long-range Coulomb interaction in 
the treatment of the few-nucleon scattering. Our further 
investigations are pointed at this purpose. 
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